Superscaling Predictions for Neutral Current Quasielastic Neutrino-Nucleus 
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The application of superscaling ideas to predict neutral-current (NC) quasielastic (QE) neutrino 
cross sections is investigated. Results obtained within the relativistic impulse approximation (RIA) 
using the same relativistic mean field potential (RMF) for both initial and final nucleons — a 
model that reproduces the experimental (e, e') scaling function — are used to illustrate the ideas 
involved. While NC reactions are not so well suited for scaling analyses, to a large extent the 
RIA-RMF predictions do exhibit superscaling. Independence of the scaled response on the nuclear 
species is very well fulfilled. The RIA-RMF NC superscaling function is in good agreement with 
the experimental (e, e') one. The idea that electroweak processes can be described with a universal 
scaling function, provided that mild restrictions on the kinematics are assumed, is shown to be valid. 
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Analyses of on-going and future experimental stud- 
ies of neutrino reactions and oscillations at intermedi- 
ate energies [l[ inevitably involve nuclear targets and 
require accurate control of nuclear effects. One way of 
taking nuclear effects into account is by directly mod- 
cling them. This approach can predict the bulk of the 
neutrino-nucleus response, but is not capable of yielding 
predictions of high enough accuracy, given the present ex- 
perimental demands. A second approach that has been 
recently proposed takes advantage of scaling ideas. In- 
deed, scaling has been extensively employed to analyze 
inclusive QE electron- nucleus scattering data 0, Q . The 
data, when appropriately organized, scale to a function 
that is not only relatively independent of the momen- 
tum transfer (scaling of the first kind) , but also indepen- 
dent of the nuclear target (scaling of the second kind). 
The simultaneous occurrence of both kinds of scaling 
is known as superscaling Q. Based on these ideas, a 
phenomenological Super Scaling Approach (SuSA) [J, Q 
can be pursued that provides a more robust way to 
inter-relate the various classes of electroweak processes 
than most direct modeling does, as long as the kine- 
matics chosen lie in the regions where scaling applies, 
i.e., QE kinematics for transferred momentum in the 
range from roughly 500 MeV/c to a few GeV/c. Within 
SuSA, one assumes that at similar kinematics both elec- 
tron and neutrino mediated inclusive scattering reactions 
share the same universal scaling function, which con- 
tains the relevant information about the initial and final 
state nuclear dynamics explored by the probe, thereby 
allowing one to provide reliable and relatively model- 
indcpcndcnt predictions for neutrino-induced processes 
employingthe (e, e') experimental scaling function as in- 
put [i 1 1, 0, @, s m to, m, m . 

To date, most applications of scaling ideas to neutrino- 
nucleus cross sections involved charged current (CC) pro- 



cesses, whose kinematics parallel the electron scattering 
case. However, the interaction of neutrinos with mat- 
ter is mediated not only by W bosons, but also by the 
neutral Z° boson. NC processes are relevant for oscil- 
lation experiments — for instance, it is expected that 
they contribute as the third most important event type 
for the MiniBooNE experiment at Fermilab [jj. As in 
the case of CC processes, predictions based on scaling 
ideas, when possible, are clearly demanded. The identi- 
fication of CC events is relatively simple via the outgo- 
ing charged lepton, similar to what happens in inclusive 
(e, e') scattering. This means that the energy and mo- 
mentum transferred at the leptonic vertex are known and 
thus the scaling analysis of CC neutrino-nucleus cross 
sections proceeds in a way identical to the electron case. 
However, in the case of NC events, the scattered neutrino 
is not detected and identification of the NC event is usu- 
ally made when i) no final charged lepton is found and ii) 
a nucleon ejected from the nucleus is detected. Even in 
the case that the nucleon energy and momentum can be 
measured, the transferred energy and momentum at the 
leptonic vertex will remain unknown. The kinematics 
of the NC process is thus different from both electron 
scattering and its CC neutrino counterpart, rendering 
the derivation of scaling less obvious. Nevertheless, the 
translation of the scaling analysis to NC processes was 
recently outlined in There it was shown that the 

superscaling analysis of NC reactions in the case of the 
Relativistic Fermi Gas (RFG) and scattering of 1 GeV 
neutrinos from 12 C is feasible. Said study showed how to 
extend the scaling analysis to NC processes. The RFG 
(e, e' ) re sponse exhibits perfect superscaling by defini- 
tion |T3 |. but it is not in accord with the magnitude or 
with the shape of the experimental scaling function. It 
has been shown that strong final-state interactions (FSI) 
are needed to describe successfully the magnitude and 
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shape of the superscaled data, introducing also small de- 
viations from the extracted superscaling behaviour. 

In this Letter, we address two crucial questions which 
arise when extending SuSA analyses to NC neutrino scat- 
tering in the QE region: i) does superscaling hold for 
NC neutrino-nucleus cross sections when strong FSI are 
present? If so, ii) can the (e, e') experimental scaling 
function be employed to predict NC cross sections, in 
spite of the intrinsic differences between the two pro- 
cesses? To answer these questions, being aware of how 
scarce NC neutrino-nucleus cross section data are, we 
use predictions from the Relativistic Impulse Approxi- 
mation (RIA) 0, 0, H GJ, [U O, based on strong rel- 
ativistic mean field potentials for both the bound and 
ejected nucleons (RIA-RMF). This model, as well as its 
corresponding semirelativistic version Q , reproduces the 
shape and magnitude of the experimental scaling curve 
extracted from QE (e, e') data, elusive for other theoreti- 
cal models. Furthermore, RIA-RMF predicts a universal 
scaling function for both electron and CC neutrino scat- 
tering [3, H, Here, we verify for the first time that 
NC QE neutrino cross sections exhibit superscaling prop- 
erties even in presence of strong FSI. Insights into the 
universal character of the scaling function, i.e., the exis- 
tence of a unique function that simultaneously describes 
QE electron, CC and NC neutrino scattering on nuclei, 
are also provided. 

In NC QE neutrino scattering an outgoing nucleon 
(mass tun) having energy En, kinetic energy Tn = 
En — to jv and angle 9k PN with respect to the momen- 
tum k of the beam is assumed to be detected. The beam 
energy e is also assumed to be known. These variables 
determine the kinematics of the process [TH, . With 
regards to the model we employ, the NC QE neutrino- 
nucleus scattering is described within the impulse ap- 
proximation (IA), where the nuclear current is written 
as a sum of single-nucleon currents. The bound nucleon 
states are given as self-consistent Dirac-Hartree solutions, 
derived within a RMF approach using a Lagrangian con- 
taining a, uj and p mesons |2pj . FSI effects are included 
by means of the same strong RMF potentials that de- 
scribe the initial bound states. A more detailed descrip- 
tion of the model can be found in [IH [ltjl [l7l l2~i| . 

As usual in scaling analyses of QE scattering, we as- 
sume the inclusive A(v,N)v'X cross section to be ob- 
tained as the integrated semi-inclusive one-nucleon (pro- 
ton or neutron) knockout A(v, v' N) X cross sections. In 
Fig. [1] we show the strong dependence of NC neutrino 
QE inclusive cross sections on the beam energy (pro- 
vided that 0k PN is fixed), and on the target selected. 
The results are obtained with the RIA-RMF model; how- 
ever a large amount of this variation is essentially due to 
the neutrino-nucleus coupling strength and the variation 
in the position of the quasielastic peak for the different 
beam energies. If superscaling holds, most of this depen- 
dence disappears when dividing these cross sections by 
the NC single-nucleon cross section given in Eq. (20) of 
[lH and plotting against the dimensionless scaling vari- 
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FIG. 1: NC QE differential cross section da / '(di?jvdf2jv) ver- 
sus the outgoing proton kinetic energy Tn for the reaction 
(y,y). The left-hand panel corresponds to 12 C at different 
incident neutrino energies e. The right-hand panel shows re- 
sults at fixed e — 1 GeV for different target nuclei. In both 
panels, the outgoing nucleon detection angle is 9k PN = 40°. 



able ip u extracted from the RFC analysis in NC kine- 
matics (see Eq. (26) in [13] for its explicit expression). 
The differences in nuclear species should also be taken 
into account by the superscaling analysis. Results for 
the so-obtained scaling function f{ip u ) are presented in 
Fig. El 
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FIG. 2: NC scaling functions corresponding to the differential 
cross sections in Fig. [1] 

In the left-hand panel of Fig. [2j one can see that first- 
kind scaling is well respected within RIA-RMF in the 
region of negative ^"-values. In other words, the large 
variations in the cross sections observed for different neu- 
trino energies, are accounted for by the single-nucleon 
part of the cross sections, which has been factored-out in 
obtaining the scaling function. Furthermore, the peak of 
the superscaling response appears approximately at the 
same point for all the kinematics. However, first-kind 
scaling is not perfect, as there is a sizeable increase in 
the height of the peaks of the curves, as well as a shift 
to ip u > for increasing beam energy. This is simi- 
lar to what is observed in RIA-RMF for the inclusive 
(e, e') case. Actually, the experimental (e, e') data do 
leave room for some breaking of first-kind scaling in the 
region of positive scaling variable. First-kind scaling is 
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very well fulfilled for electron, CC and NC cases in the ab- 
sence of FSI @, H d GJ, HI. Therefore, the breakdown 
of scaling in Fig. 2 must be ascribed (within IA) to FSI. 
In the plane- wave limit, the dependence of the cross sec- 
tion on the energy of the outgoing nucleon comes mainly 
from kinematical effects that are taken into account in 
the scaling analysis. However, FSI involve a redistribu- 
tion of strength that depends on the energy of the final 
nucleon. In other words, FSI introduce an additional, 
non-kinematical, dependence of the cross section on Tpj. 
If the kinematics of the process are such that the range 
of energies of the ejected nucleon depends strongly on 
the beam energy, the nucleon will be subject to different 
FSI for each e, and a visible breakdown of first-kind scal- 
ing will show up. This is what happens for 9kp N = 40°, 
where there is a strong shift of the position of the peak of 
the cross section with incoming beam energy. However, 
for those kinematics for which the range of Tjy remains 
approximately the same when considering different beam 
energies, first-kind scaling is obtained even with FSI in- 
cluded, as FSI effects on the knockout nucleon are similar 
for different beam energies. 





1800 




1600 




1400 




1200 


> 






1000 


5 




z 


800 


1- 


600 




400 




200 








1 GeV 

1.5 GeV 

2 GeV 




FIG. 3: Relationship between 6k PN and Tjv for the NC (y,p) 
reaction on free protons at rest. The different curves show 
the results for three beam energies, 1, 1.5 and 2 GeV. 

Incidentally, in Fig. 2 we also observe that f(ip u ) is, 
for pure kinematical reasons, a bivalued function of the 
scaling variable ip u , as the same value of ip u is obtained, 
at fixed beam energy and nucleon angle, for two different 
values of the o utg oing nucleon energy. In the absence of 
FSI (as in ref. [13|). superscaling is a good approximation 
and the two values of the superscaling function for these 
ijj u are nearly equal. When FSI are present, and if the 
kinematics prevents superscaling, the bivalued nature of 
the superscaled function is revealed. 

In order to understand for what kinematics good scal- 
ing of the first kind is reached even in presence of strong 
FSI, we look at the case of free nucleons. Fig. [3] shows 
how 9kp N and Tjv are related due to energy and momen- 
tum conservation for several beam energies. For bound 
nucleons, neglecting Fermi motion, the cross section will 
be peaked at approximately the same Tjv value. From the 
figure one sees that the range of Tjv spanned at fixed 9k PN 
for varying e is reduced for large angles and thus scaling 



of the first-kind will be much better obeyed. In general, 
smaller angles show larger first-kind scaling violations, 
while larger angles exhibit almost perfect first-kind scal- 
ing [22| . Note that this result comes through purely kine- 
matical reasoning and thus is model-independent to the 
extent that the cross section can be described within IA. 

Results for scaling of the second kind are presented in 
the right-hand panel of Fig. [2] The superscaling func- 
tions obtained for several nuclei are almost identical, in 
spite of the strong difference in magnitude of the cor- 
responding cross sections (cf. Fig. [I]). That is, the de- 
pendence on the nuclear species is well accounted for by 
the superscaling analysis. Scaling of second kind is seen 
to be very robust, thereby opening up a means of taking 
into account nuclear effects for different nuclei employing 
superscaling ideas. 

The superscaling properties exhibited by NC QE 
neutrino-nucleus scattering suggest exploring the validity 
of the universal character of the scaling function for in- 
clusive electroweak processes on nuclei, using either elec- 
trons or CC and NC neutrino probes. To the extent that 
this universality holds, the phcnomcnological SuSA ap- 
proach formerly applied to predict CC neutrino-nucleus 
cross sections could also provide reliable, largely model- 
independent predictions for NC processes. 




FIG. 4: NC scaling function evaluated within the RIA-RMF 
approach for 40 and 60 degrees at two beam energies, com- 
pared with the averaged experimental function, together with 
a phenomenological parameterization of the data (solid). 

In order to study whether or not this universality as- 
sumption holds also for NC processes, in Fig. 0] we com- 
pare the RIA-RMF NC superscaling function with the 
averaged QE experimental function obtained from the 
analysis of (e, e') data, together with a phenomenologi- 
cal parameterization [g, H [23| ■ The RIA-RMF superscal- 
ing function has been plotted for two values of 9kp N for 
which scaling of first-kind is well fulfilled (60°) or not- 
so-well (40°). Results are shown for two beam energies. 
As seen, the model gives rise to a NC scaling function 
that follows closely the behaviour of the (e, e') function 
and the bivalued behaviour of the superscaling function 
is hardly visible whenever superscaling is well respected 
(60°). For the case of 40°, breakdown of first-kind scal- 
ing is clear, the departure from the SuSA (e, e') response 
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visible and the bivalued nature of the NC superscaled 
function is enhanced. We notice that all curves would co- 
incide if superscaling was exactly fulfilled in both NC and 
(e, e') cases. Since the (e, e') and NC scaling curves are 
obtained under rather different kinematical situations, 
the scaling curves depart from each other when super- 
scaling is not a good approximation. This supports the 
assumption that, under proper kinematics restrictions, a 
universal QE scaling function exists which is valid not 
only for inclusive electron and CC neutrino reactions, as 
seen in @, [1], but also for NC processes. 

In summary, we have established sufficient conditions 
under which a universal superscaling function could be 
applied both to electron and neutrino (CC or NC) inclu- 
sive scattering. These conditions refer to the fact that 
the kinematics must be such that the range of energies 
spanned by the ejected nucleon is nearly independent 
of the incoming neutrino energy. This happens, for in- 
stance, when the angle of the ejected nucleon with re- 
gards to the beam is larger than roughly 50°, which hap- 
pens to be the region where the cross section integrated 
over angles has larger values. In such a case, first-kind 
scaling is well respected at the 10% level even in the pres- 
ence of strong FSI, and the good comparison with the 
experimental (e, e') scaling function gives us confidence 
that SuSA can be extended to predict NC QE neutrino 



cross sections. 

We also note that even though we have illustrated this 
study within the RIA-RMF model that contains strong 
FSI and is quite successful in reproducing the experimen- 
tal electron scattering scaling function, the kinematical 
conditions that grant the validity of SuSA are model in- 
dependent when the IA can be safely applied, that is 
under QE kinematics with neutrino beam energies from 
- 500 MeV up to a few GeV. 
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